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^  Abstract 


An  abstract  approximation  framework  and  convergence  theory  for  the  identification  of 
thermoelastic  systems  is  developed.  Starting  from  an  abstract  operator  formulation  consisting 
of  a  coupled  second  order  hyperbolic  equation  of  elasticity  and  first  order  parabolic  equation  for 
heat  conduction,  well-posedness  is  established  using  linear  semigroup  theory  in  Hilbert  space,  and  a 
class  of  parameter  estimation  problems  is  then  defined  involving  mild  solutions.  The  approximation 
framework  is  based  upon  generic  Galerkin  approximation  of  the  mild  solutions,  and  convergence 
of  solutions  of  the  resulting  sequence  of  approximating  finite  dimensional  parameter  identification 
problems  to  a  solution  of  the  original  infinite  dimensional  inverse  problem  is  established  using 
approximation  results  for  operator  semigroups.  An  example  involving  the  basic  equations  of  one 
dimensional  linear  thermoelasticity  and  a  linear  spline  based  scheme  is  discussed  and  numerical 
results  indicating  how  our  approach  might  be  used  in  a  study  of  damping  mechanisms  in  flexible 
structures  are  presented.  /'  1 
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1.  Introduction 


In  this  paper  we  develop  an  abstract  approximation  framework  and  convergence  theory  for  the 
identification  of  abstract  linear  thermoelastic  systems.  The  inclusion  of  thermal  effects  in  the 
dynamics  of  flexible  structures  has  recently  received  an  increased  amount  of  attention  as  an  area  of 
research  as  a  result  of  the  problem  of  solar  heating  of  large  flexible  spacecraft  and  the  effect  that  it 
can  have  on  the  structure’s  vibrational  modes  and  the  subsequent  design  of  efficient  and  effective 
control  laws. 

The  approach  we  take  here  is  somewhat  different  than  the  traditional  one  usually  taken  by  other 
authors  working  in  the  area  of  thermoelasticity  (see,  for  example,  [4],  [5],  [9],  [13]).  Rather  than 
starting  with  the  mathematical  formuation  of  the  basic  laws  of  continuum  mechanics  and  thermo¬ 
dynamics  and  then  linearizing  to  arrive  at  the  so  called  basic  equations  of  linear  thermoelasticity, 
we  consider  an  abstract  operator  formuation  of  these  basic  equations  set  in  appropriately  cho¬ 
sen  infinite  dimensional  spaces.  More  precisely,  we  consider  an  abstract  second  order  equation  of 
elasticity  coupled  with  a  first  order  abstract  parabolic  equation  to  describe  the  dynamics  of  heat 
conduction.  We  then  rewrite  this  system  as  an  equivalent  abstract  first  order  system  set  in  an 
appropriately  constructed  product  Hilbert  space.  We  establish  well-posedness  using  results  from 
linear  semigroup  theory  (i.e.  the  Lumer- Phillips  Theorem)  and  then  proceed  to  define  a  class  of 
inverse  or  parameter  identification  or  estimation  problems. 

Our  approximation  results  are  in  the  spirit  of  the  treatment  by  Banks  and  Ito  in  [1].  We  develop 
an  abstract  approximation  framework  via  generic  Galerkin  approximations  to  mild  solutions  of  the 
abstract  state  equations.  A  sequence  of  finite  dimensional  approximating  parameter  identification 
problems  result.  Using  an  approximation  and  convergence  result  from  linear  semigroup  theory  (the 
Trotter-Kato  Theorem;  stability  and  resolvent  convergence  imply  semigroup  convergence)  we  are 
able  to  establish  that  under  reasonable  and  readily  verifiable  assumptions,  solutions  to  the  finite 


dimensional  identification  problems,  in  some  sense,  approximate  a  solution  to  the  original  infinite 
dimensional  inverse  problem. 

We  note  that  in  a  certain  sense,  the  results  that  we  present  here  subsume  the  framework  developed 
in  [1].  Indeed,  Banks  and  Ito  first  treat  first  order  abstract  parabolic  systems  and  then  consider 
second  order  systems.  The  abstract  thermoelestic  systems  of  interest  to  us  here  consist  of  one 
equation  of  each  type  appropriately  coupled.  By  zeroing  out  the  couDli"?  ( which  the  framework 
we  develop  below  can  handle)  in  our  treatment,  we  obtain  most  of  the  results  for  first  and  second 
order  systems  derived  separately  by  Banks  and  Ito  in  [1],  simultaneously. 

In  order  to  demonstrate  that  in  addition  to  being  theoretically  sound,  our  approximation  frame¬ 
work  is  useful,  numerically  and  computationally  feasible,  and  performs  well  in  practice  (or  at  least, 
simulation)  we  apply  our  general  approach  to  the  basic  equations  of  linear  one  dimensional  ther¬ 
moelasticity.  This  system  consists  of  a  one  dimensional  wave  equation  describing  the  longitudinal 
or  axial  vibrations  of  a  slender  rod,  coupled  with  the  heat  equation  together  with  appropriate 
boundary  and  initial  conditions.  We  illustrate  how  a  spline  based  scheme  can  be  constructed  so 
as  to  conform  to  our  framework  and  indicate  how  it  might  be  used  in  a  study  of  dissipation  or 
damping  mechanisms  in  a  flexible  structure.  Numerical  results  carried  out  on  a  Cray  XMP/48 
supercomputer  are  presented  and  discussed. 

An  outline  of  the  remainder  of  the  paper  is  as  follows.  In  the  second  section  we  define  the  abstract 
thermoelastic  system,  establish  its  well-posedness,  and  state  precisely  the  class  of  identification 
problems  which  are  of  interest  to  us.  The  approximation  framework  and  convergence  theory  are 
developed  in  the  third  section,  and  the  example  and  numerical  studies  are  discussed  in  Section  4. 
Section  5  contains  a  brief  summary  of  our  findings  and  some  concluding  remarks. 
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2.  Abstract  Linear  Thermoelastic  Systems 

In  this  section  we  define  an  abstract  linear  thermoelastic  system,  and  demonstrate  the  existence, 
uniqueness,  and  continuous  dependence  of  solutions  via  linear  semigroup  theory.  We  then  define 
an  associated  identification  or  inverse  problem. 

Let  Q  be  a  metric  space  and  let  Q  be  a  compact  subset  of  Q.  The  set  Q  will  be  known  as  the 
admissible  parameter  set.  For  j  =  1,2  let  {H  {■,  •)■  be  real  Hilbert  spaces  and  let  {V),||’||j} 
be  reflexive  Banach  spaces.  We  assume  that  for  j  =  1,2,  Vj  is  densely  and  continuously  embedded 
in  Hj  with  \<p\-  <  Hj  ||<^|| .  ,y?  G  Vj.  We  let  V*  denote  the  continuous  dual  of  Vj.  Then  with  Hj  as 
the  pivot  space,  we  have  Vj  <-*  Hj  =  HJ  |— ►  V*,  with  Hj  densely  and  continuously  embedded  in 
V*.  We  denote  the  usual  operator  norm  on  V*  by  || '|| , J  =  1,2.  In  the  usual  manner,  {•,  -)j  is 
understood  to  denote  both  the  inner  product  on  Hj  and  the  duality  pairing  on  V*  X  Vj,  for  j  —  1,2. 

For  each  ?  G  Q  we  consider  what  we  shall  refer  to  as  an  abstract  linear  theormoelastic  system 
given  by 

u  (t)  +  C  (q)u{t)  +  Ai(q)u(t)  +  L{q)’ 6(t)  =  f  (t\q) ,  t>  0  (2.1) 

6(t)  +  A2  (q)O(t)  -  L(q)u(t)  =  g  {t',q),  t  >  0,  (2.2) 

u  (0)  =  «o  (?) ,  u  (0)  =  n0  (?) ,  0  (0)  =  60  (?) ,  (2.3) 

where  for  each  t  >  0,u(t)  G  H\  and  6  (t)  G  fh-  We  assume  that  for  each  q  G  Q,C  {q),Ai  (?)  G 
C(Vl,Vl’),A2(q)eC(V2,V{),  L(q)e  C(VuV{),uo  (?)  €  Vj,  v0  (?)  G  H\,0o  (?)  €  H2,  and 
/(•;?)  G  Li  (0 ,T;Hi),  ?(•;?)  €  L\  (0 ,T;  H2),  for  some  T  >  0.  The  operator  T(?)*  G  L(V 2,V{)  is 
defined  by 

(L{q)*  xl>,y)x  =  (L(q)q?,ip)2  ,ip  G  V2,<p  G  Vv  (2.4) 

We  shall  also  require  the  following  further  assumptions. 

(A)  (Symmetry)  For  each  ?  G  Q  the  operator  A\(q)  is  symmetric  in  the  sense  that 


(Ai  (?)??,  V’)i  =  (^1  (?)i/>,?>}i  ,¥>,  ^  G  Vi. 


(B)  (Uniform  Coercivity)  For  j  —  1,2  there  exist  constants  atj,(3j  G  R,  independent  of  q  G  Q,  with 
/ 3j  >  0  such  that 

Mi  (</)<?,  ¥>).,•  +  oy  |v>|j  >  fa  \\(p\\2.  ,q>  G  V-. 

(Note  that  without  loss  of  generality  we  may  assume  that  Oj  >  0,j  =  1,2.) 

(C)  (Uniform  Boundedness)  There  exist  positive  constants  7j,ji  =  1, 2,  <5,  and  p,  independent  of  q  G  Q, 
for  which 

<  7j  €  Vj,j  =  1,2 

!IM9MI2»  <  pIMIi  ,9  e  Vi. 

(Note  that  the  final  bound  above  implies  that  ||T  (</)*  <  p  ||0||2  ,  ip  £  U2,  as  well.) 

(D)  (Nonnegativity)  For  each  q  £  Q  the  operator  C\q)  is  nonnegative  in  the  sense  that  (C  ( q )  > 

0,  v3  £  Vj. 

(E)  (Continuity)  For  G  Ui  and  ip  £  V 2  the  mappings  9  -*  Ai  {q)<p,q  — *•  C  (<7)9,9  — *•  L(q)*  ip,  q  — <• 
A2  (9)^,  and  9  — >  L(q)(p  are  continuous  from  Q  C  Q  into  Vj*  or  V2*  (which  ever  is  appropriate). 
Also,  the  mappings  9  — ►  uo(q),q  —>  v0(q),  and  q  — ►  60  {q)  are  continuous  from  Q  C  Q  into 
Vi,Hi,  and  /f2,  respectively,  as  are  the  mappings  q  — *•  f(t;q )  and  q  — >  g{t\q)  into  7/j  and  772, 
respectively,  for  almost  every  t  G  [0,T]. 

(Note  that  the  two  “strong”  continuity  assumptions  on  the  operators  L(q)  and  L(q)m  can  be 
replaced  by  the  single,  somewhat  stronger  continuity  assumption  that  the  mapping  q  —*  L{q)  is 
continuous  from  Q  C  Q  into  £(V i,  V2*)  -  i.e.  q  — *  L(q)  is  continuous  with  respect  to  the  uniform 
operator  topology  on  C{V\ ,  V2*).) 

(F)  (Uniform  Domination)  There  exist  functions  /o,  go  E  Li(0,T),  independent  of  q  G  Q,  for  which 
1/ (t',q)\i  <  fo  (t)  and  \g{t;q)\2  <  g0  (<),  a.e.  t  G  [0,T]  and  every  q  G  Q. 


Remark  A  comment  regarding  our  referring  to  the  system  (2.1)  -  (2.3)  as  an  abstract 
thermoelastic  system  would  be  in  order  here.  If  one  thinks  of  (2.1)  as  an  abstract  elasticity  equation 
and  (2.2)  as  an  abstract  heat  equation,  then  by  making  some  rather  simple  assumptions  concerning 
thermal  stress,  and  applying  some  basic  mechanical  and  conservation  of  energy  principles  from 
continuum  and  thermal  mechanics,  one  could  justify  the  form  of  the  coupling  between  the  two 
equations  (i.e.  the  operators  L(q)  and  L(q)m  ).  Typically  however,  the  basic  equations  (both  linear 
and  nonlinear)  of  thermoelasticity  are  formulated  from  first  principles  using  nontrivial  mathematical 
formulations  of  the  basic  laws  of  mechanics  and  thermodyamics  (see,  for  example,  [4],  [9],  or  [13]). 
In  light  of  this,  we  justify  our  reference  to  (2.1)  -  (2.3)  as  an  abstract  thermoelastic  system  by 
simply  stating  that  the  thermo-mechancial  systems  of  interest  to  us  here,  and  in  particular  the 
basic  equations  of  linear  thermoelasticity  in  one  and  higher  dimensions  (see  [4], [5],  and  [16]),  can, 
via  appropriate  identifications,  be  put  in  this  abstract  form. 


We  demonstrate  the  well-posedness  of  the  system  (2.1)  -  (2.3)  for  each  q  £  Q  by  first  rewriting 
it  as  an  equivalent  first  order  system  in  an  appropriate  product  Hilbert  space  and  then  applying 
results  from  linear  semigroup  theory. 

Let  X  be  the  Banach  space  defined  by  A'  =  V\  X  H\  X  Hi  with  norm  |  •  |x  given  by  |(y?,  ip ,  77) |  v  = 
(iMIi  +  {i’ll  +  M2)  •  For  each  q  E  Q  let  X(q)  denote  the  Hilbert  space  which  is  set  equivalent 

to  X  and  which  is  endowed  with  the  inner  product  (-,-)q  given  by 


(Ml,  01  =  (^1  («)V>1  ,<P2)i 


+«i  {'Pu&h  +  (01  ,fa)i  +  (91,92)2 , 

for  E  X,  i  =  1,2.  We  denote  the  norm  on  X(q)  induced  by  the  inner  product  {•,•)  by 

|  j  and  note  that  it  is  immediately  clear  that  assumptions  (B)  and  (C)  imply  that  the  norms  |  •  |_\- 
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and  |  •  |9  are  equivalent,  uniformly  in  g  for  q  £  Q.  That  is,  there  exist  positive  constants  m  and  M, 
independent  of  q  £  Q,  for  which 


m  H9  <  l-lx  £  M  I 


(2.5) 


For  each  q  £  Q  define  the  operator  4(g)  :  Dom(A(q))  C  X  (q)  —>  X  (q)  by 


A(q)  = 


0/0 

-Ai(q)  —C  (q)  —  L  (q) 

0  L(q)  -42(g)J 


(2.6) 


Dom(A  (q))  =  €  Vx  x  Vt  x  V2  :  4(g)(y>,  ^,9)  £  X)  .  (2.7) 

Theorem  2.1  For  each  q  £  Q  the  operator  4(g)  defined  in  (2.6),  (2.7)  above  is  the  infinitesimal 
generator  of  a  Co-semigroup,  {5  (t\q)  :  t  >  0}  of  bounded  linear  opearors  on  X(q)  (and  therefore 
X  as  well  ). 

Proof  The  result  will  follow  from,  for  example,  Showalter  [12],  Theorem  IV.4.C,  if  we  can  show 
that  for  some  w  €  R  the  operator  A  ( q )  -  ujI  is  dissipative  and  the  operator  A  -  A  ( q )  is  surjective 
for  some  A  >  u>.  Toward  this  end,  for  x  =  £  Dom(A(q))  we  have 

(A(q)x,x)g  =  ((0,  —  Ai  (g)y>-C(g)0-  L(q)'  q,L{q)xp  -  A2  (q)  tj),  (<p,  0,  Tj))g 
=  (9)MMi  +  tti  (V»,V?)i  -  (Ai  {q)<p,tj>) j  -  (C  (g)V>», 

~  (£(?)’ 7?,^),'+  (L  (q)ip,  t/)2  -  (A2(q)ri,il)2 
=  «i  (</>,  V»)i  “  (C  (9)  Vb  Mi  “  (^2  (9)  9, 9)2 

<  yMi  +  YH?-/32||9ll?+o,|9g 

<  yy  IMIi  +  Y  Mi  +  «2 19I2 

-  {<Al  (9)^, ^>1  +«i  Mi}  +  |lMi  +  <*2  M2 
<w|x|J, 

where  u>  =  max  j  ,  ^jS«2  }•  Recall  that  we  assumed,  without  loss  of  generality,  that  01,02  >  0, 
and  note  that  when  ai  =  a2  =  0,  we  have  u  =  0  and  that  the  operator  4(g)  is  dissipative. 
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In  order  to  eliminate  some  rather  tedious  technical  details,  we  assume  for  the  remainder  of  the 


proof  that  ai  =  o2  =  0  (and  therefore  u>  =  0).  The  arguments  to  follow  remain,  in  substance, 
essentially  unchanged  in  the  more  general  case  of  0:1,02  >  0.  Let  A  >  0,  let  (0,0,7?)  G  X  and 


consider 


or  equivalently 


(A  -  A(q))(u,v,6)  =  (0,0,7?); 


Xu  —  v  =  0 


A  v  +  Ax  (q)  u  +  C  (q)  v  +  L  (q)’  0  =  0 


A  6  -  L{q)v  +  A  2  (q)  0  =  7?. 


(2.10) 


Solving  (2.8)  for  v  and  then  substituting  into  (2.9)  and  (2.10)  we  obtain 


(A2  +  A  C  ( q )  +  Ax  («?))  u+L  (q)m  9  =  0  +  (C  (?)  +  A)  0 


(2.11) 


-L(q)u  +  (l  +  A  xA2{q))0  =  X  1  (7?+  L(q)ip). 


(2.12) 


Let  H  =  H\  X  H 2  be  endowed  with  the  inner  product  (■,■)  given  by 


($,*)  =  +  (*2, *2)2, 


(2.13) 


for  $  =  ($i,$2),4'  =  (^1,^2)  G  H,  and  let  V  =  Vj  x  V2  be  endowed  with  the  norm  ||-||  given  by 


11*11  =  (ll*.|li  +  IM2)''J, 


(2.14) 


for  $  =  ($1,^2)  G  V.  Then  H  is  a  Hilbert  space,  V  is  a  reflexive  Banach  space  and  V  <— ►  H  ' >  V*, 
with  the  embeddings  dense  and  continuous.  Define  the  operator  A\  ( q )  €  C  (V,  0*)  by 


[  A2  +  AC  (<?)!  .!,(<?)  L(q)'  1 

Al9)_[  -L(q)  (1+A-M2(9))J  ' 


(2.15) 


Then  for  <J>  =  ( <1* i ,  *I>2 )  £  V,  it  is  easily  shown  that 

(A\  (</)  4>,  <I>)  >  J,  \\<t>l  ||'{  +  A-1  *||*2||:j 


(2.10) 


11*11  , 

where  f3\  =  min  {/Jj, A_1/3o}  >  0.  It  therefore  follows  that  A\(q)  is  an  isomorphism  from  \'  onto 
Vm  (see,  for  example,  Tanabe  [15],  Theorem  2.2.2).  Recalling  (2.8),  (2.11)  and  (2.12)  and  setting 


(u,8)  —  A\(q)  (v  +  [C  (q)  +  A)  *p,  A  1  (i)  +  L  (q)  ^)) 


(2.17) 


and 

v  =  A  u  —  <p,  (2.18) 

we  obtain  that  A  -  A(q)  is  a  surjection,  and  the  theorem  is  proved. 

For  each  q  €  Q  define  x  (•;</)  6  C  ([0. 7'] :  .V  )  by 

t 

x(t-.q)  =  -S'  (/;</)  x,)  (q)  +  J  S(l  -  *\q)  F(.s:q)<Ls.  (2.19) 

o 

for  t  €  [0 ,T]  wliere  {S(t\q)  :  I  >  0}  is  the  semigroup  generated  by  A  (q)  ,xo(q)  =  (  uq  (</) .  t'o  (q)  ■ 
#o  (q))  €  A',  and  F{t\q)  =  (0,f(t;q).g(t;q)),  a.e.  t  €  [0,  T].  We  shall  refer  to  the  A’ -valued 
function  x(-\q)  =  (u  (•;  q) ,  v  (•;  q)  ,9  (•;  q))  given  by  (2.19)  above  as  the  unique  mild  solution  to  the 
abstract  theormoelastic  system  (2.1)  -  (2.3). 

We  now  proceed  to  define  an  identification  problem  corresponding  to  (2.1 )  -  (’>  3).  Let  Z  denote 
an  observation  space.  For  i  =  1,2, . . .  , v  and  z  €  Z,  let  4>,  (•,•;£)  :  X  xQ  —  R+  denote  a  continuous 
map  from  X  X  Q  C  A'  X  Q  into  the  nonnegative  real  numbers.  We  consider  the  following  parameter 
estimation  or  inverse  problem. 


(ID)  Given  observations  {r,-}*! j  C  X  Z  at  times  { 2 , } j  C  x  [0.7’]  determine  parameters  q  £  Q 

i=i  i=t 

which  minimize 

V 

J(q)  =  (220) 

t=i 
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where  for  each  q  £  Q  and  t,  £  [0,  T],  x  (/,;  q)  is  given  by  (2.19). 

We  do  not  address  the  question  of  t ho  existence  of  a  solution  q  to  problem  (ID)  here,  but  rather 
defer  this  discussion  until  the  next  section.  In  the  section  to  follow,  the  existence  of  a  minimizer 
for  the  functional  J  given  by  (2.20)  is  established  as  a  consequence  of  our  approximation  theory. 


3.  Approximation  and  Convergence 

Using  a  standard  Galerkin  technique  we  construct  a  sequence  of  finite  dimensional  approximations 
to  the  abstract  thermoelastic  system  (2.1)  -  (2.3).  A  sequence  of  approximating  identification 
problems  result.  Using  a  well  known  approximation  theoretic  result  for  linear  semigroups  (i.e.  the 
Trotter-Kato  theorem)  we  establish  that  a  solution  to  the  infinite  dimensional  parameter  estimation 
problem,  (ID),  exists,  and  moreover,  that  it  is  in  some  sense  approximated  by  solutions  to  the  finite 
dimensional  approximating  identification  problems. 

For  j  =  1,2  and  for  each  n j  =  1.2....  let  //"’  be  a  finite  dimensional  subspace  of  II :  with 
//"’  C  Vj ,  for  all  Uj.  Let  /J"’  :  II }  —  //"'  denote  the  orthogonal  projection  of  II j  onto  //"’ 
computed  with  respect  to  the  (•,•)  ■  inner  product.  We  shall  require  the  following  approximation 
assumption. 

(G)  (Approximation)  For  j  =  1,2 


lim 


v~e  Vj. 


Note  that  the  dense  and  continuous  embedding  of  Vj  in  II }  together  with  assumption  (G)  also  yield 
that  lim  IP"' q>  -  q>\  =  0,  q>  £  //,,  j  =  1, 2. 

n,  — »oo  *  1 3 

For  each  q  £  Q  we  define  the  operators  A"'  (q)  £  C  (//"’)  ,J  =  1,2,  C1  (q)  £  C  (//"’).  and 
Ln  (q)  £  C(H^' ,  //2"3 ),  n  =  ( 7i i ,  n2 ),  using  standard  Galerkin  approximation.  More  precisely,  for 
j  =  1,2,  and  q>n'  £  //"',  we  set  A™’  (q)q>n>  =  4'n’  £  //"'  where  is  defined  via  the  Riesz 
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Representation  Theorem  on  t he  Hilbert  space  //”’  as  (Aj{q)^pn’,\n’)]  =  {vn> , \”'  )j ,  G  H  . 
Similarly  we  set  Cni{q)<fni  =  ipn>  where  ipni  £  //”'  satisfies  (C  (?)  i^71' , \”!  )t  =  (t’"1 , \n‘)i  ■ 

Xni  G  and  £n(?)y>n 1  =  t/V12  where  G  7722  satisfies  (£  (?)  <^n> ,  xn=)2  =  (xpn\ \"2)2, 

\n2  G  7/j12.  We  define  the  operator  Ln  (q)'  £  C  (T/"2,  Z/"1 )  to  be  the  Hilbert  space  adjoint  of  the 
operator  £n  (?). 

VVe  set  Uq‘  (?)  =  P"1  u0  (?) ,  «£ 1  (?)  =  P\'VQ  (?)  and  ^q2  (?)  =  (?),  and  set  /">  (t;g)  = 

P”'/(t;?)  and  gn-{t\q)  =  P22g(t;q)  for  almost  every  t  G  [0,7’].  VVe  then  consider  the  finite 

dimensional  system  of  ordinary  differential  equations  in  7/n  =  //”’  x  77."2  given  by 

iin  (0  +  Cn'  (?)  (t )  +  4[“  (?)un  (0  +  £"  (?)*  0"  (t) 

(3.1) 

=  /n'(h?),  <>0, 

0n  (<)  +  4?*  (?)  0"  (t)  -  Ln  (?)  ti"  (0  =  ?n2  (f,q) ,  *  >  0,  (3.2) 

(0)  =  u"‘  (?) ,  iin  (0)  =  c0n’  (?) ,  0"  (0)  =  C  (?)  •  (3.3) 

VVe  next  proceed  to  rewrite  (3.1)  -  (3.3)  as  an  equivalent  first  order  system.  For  each  ni,n2  = 

1,2,...  and  n  =  (ni,n2)  let  A'n  =  //[“  X  II x  II27  be  considered  as  a  subspace  of  the  Banach 
space  X,  and  for  each  q  £  Q  let  A’"  (?)  =  Xn  be  considered  as  a  subspace  of  the  Hilbert  space 

X(q).  For  each  ?  £  Q  let  An  (?)  G  C  (A'n)  be  given  by 

0  /  O' 

An  (q)  =  -A?'(q)  -Cn>  (?)  -£"(?)*  , 

0  Ln  (?)  -A?2(?)_ 

let  Fn{V,q)  =  (0 ,/n‘  (/;?), g7*2  (*;?)),  a.e.  I  £  [0,7],  and  let  x"(?)  =  ( 1  (?),!>£'  (?Mo2  (?))• 
Setting  xn  (t)  =  (un  (<) ,  tin  (<) , 07'  (< )),  we  rewrite  (3.1)  -  (3.3)  as 

in  (0  =  A"  (?)XB(0  +  FB  (<;?),  *>0  (3.4) 


xn  (0)  =  x£  (?) .  (3.5) 

The  solution  to  the  initial  value  problem n  (3.4),  (3.5)  is  given  by 

t 

xn  (t;q)  =  Sn(t\q)x%  (?)  +  j  Sn  (t - s\q)  Fn  (s;?) ds  (3.6) 

o 
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for  t  £  [0,  7] ,  where  the  semigroup  of  bounded  linear  operators  on  Xn  (or  Xn{q) ),  {S’71  (t;  q)  :  t  >  0}, 
is  given  by 


Sn{t;q)  =  exp(tAn(q)),  t>  0. 

We  define  a  sequence  of  finite  dimensional  approximating  identification  problems  as  follows. 

1/  1/ 

(ID71)  Given  observations  {z,}"=1  C  X  Z  at  times  {ii}^=1  C  X  [0, T],  determine  parameters  qn  6  Q 

1=1  i=i 

which  minimize 

U 

Jn(q)  =  ^$i(arn(<,-;g),9;*«) 

i  =  l 

where  for  each  q  £  Q  and  U  £  [0,7]  ,  xn  (f,;<?)  is  given  by  (3.6). 

The  type  or  kind  of  convergence  results  that  we  are  about  to  summarize  and  discuss  here  have 
been  presented  in  detail  in  a  number  of  places  in  the  literature  in  the  context  of  a  variety  of 
identification  problems  for  various  types  of  distributed  parameter  systems  (see,  for  example  [1] 
and  [2]).  In  the  context  of  the  problems  (ID)  and  (ID71)  defined  above,  our  general  convergence 
framework  and  theory  takes  the  following  form. 

Using  standard  continuous  dependence  results  for  ordinary  differential  equations  (see,  for  exam¬ 
ple,  [8],  Theorem  1.3.2)  it  is  easily  argued  that  the  map  q  — >  xn(t\q)  is  continuous  from  Q  C  Q 
into  Xn  for  each  t  £  [0,7],  and  each  n  £  Z+  x  Z  +  .  It  follows  therefore,  that  the  map  q  — >  Jn(q )  is 
continuous  for  each  n.  Consequently,  since  Q  has  been  assumed  to  be  a  compact  subset  of  <2,  the 
existence  of  a  solution  qn  to  problem  (IDn)  for  each  n  is  assured. 

Form  a  directed  set  from  Z+  X  Z+  in  the  canonical  way.  That  is,  for  rn  =  (mi, m2),  n  = 
( 77 1 ,  ri2 )  €  Z+  X  Z+,  we  say  m  <  n  if  and  only  if  mi  <  rii  and  m2  <  112-  Consider  next  the  following 
proposition. 
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Proposition  3.1  If  {p"}  is  a  net  from  the  directed  set  Z+  X  Z+  into  Q  C  Q  with  p"  — >  p,  then 
xn(t;qn)  — ►  x(t-,q)  for  each  t  £  [0,T]  where  xn  and  x  are  given  by  (3.6)  and  (2.19),  respectively. 


It  is  straight  forward  to  argue  (see  [2],  Chapter  III,  Section  1)  that  if  Proposition  3.1  can  be 
verified,  then  the  fact  that  Q  is  a  compact  subset  of  Q  implies  that  the  net  {p"}  formed  from  the 
solutions  to  the  problems  (ID"),  admits  a  convergent  subnet,  {pm},  whose  limit,  {</},  is  a  solution 
to  Problem  (ID).  Thus  if  Proposition  3.1  can  be  verified,  we  have  established  the  existence  of  a 
solution  to  problem  (ID),  and  the  fact  that  it  is  in  some  sense  approximated  by  the  solutions  to 
the  finite  dimensional  identification  problems  (ID"). 

Since  Galerkin  approximation  was  used  to  construct  the  operators  An-'  (q)  ,j  =  1,2,C"'  (p),  and 
Ln  (p),  the  same  arguments  used  to  establish  that  the  operator  A  (q)-ul  is  dissipative  in  the  proof 
of  Theorem  2.1,  yield  that  the  operators  An  (p)  —u>I  are  dissipative,  for  all  n  6  Z+  X  Z+  and  q  G  Q 

with  u  independent  of  n  and  q.  It  follows,  therefore,  that  the  semigroups  {Sn(t;q)  :  t  >  0}  a  e 

uniformly  exponentially  bounded,  uniformly  in  n  and  q.  That  is, 

\Sn(V,q)\q  <e“\  t>  0, 

or,  recalling  (2.5),  that 

M 

\Sn{f,q)\x<—e“l,  l>0. 

m 

Thus,  using  assumptions  (E)  -  (G),  Proposition  3.1  will  follow  immediately  if  we  can  show  that 

lim  \Sn(t;qn)Pnx-  S(f,q)x\x  =0  (3.7) 

n—*oo 

for  each  x  £  X  and  t  £  [0,T]  whenever  qn  — ►  p,  where  the  projection  like  mapping  Pn  :  X  —> >  Xn 

is  given  by  Pn  (y>,V>,??)  =  (P"V,  P^V)-  for  (<^,Vv/)  €  X. 
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We  establish  (3.7)  using  a  version  of  the  well-known  Trotter-Kato  semigroup  approximation 
theorem  given  in  [2]  (Theorem  II. 1.14  on  page  40).  According  to  this  theorem,  the  convergence  in 
(3.7)  will  follow  once  we  have  shown  that  for  some  A  >  u> 


lim 

n— ►  oo 


(A  -  A"  (<7n))-1  Pnx  -  Pn{\-  A  (q)) 


-l 


=  0 


(3.8) 


for  each  x  £  X  wherever  qn  — »  q. 

In  light  of  (2.5)  and  assumption  (G),  the  convergence  stated  in  (3.8)  will  follow  if  we  can  argue 
that 


lim 

71— *00 


(A  -  An  (qn))~l  Pnx  -  (A  -  A{q))~l  x 


-l 


=  0 


(3.9) 


for  each  x  £  X  and  some  A  >  u>,  wherever  qn  — *  q.  For  simplicity,  we  again  assume  that  =  a2  =  0 
and  therefore  that  u>  =  0.  The  proof  in  the  more  general  case  when  c*i  =  o2  >  0,  is  essentially  the 


same. 


Suppose  that 


(A-  A(q))(u,n,0)  =  (V’.’M), 


(3.10) 


and 


(3.11) 


(A  -  An  (</n))(«n',rni,r)  =  Pn(*>,<M) 

=  (prv,pri>,p?Ti)- 

We  shall  show  that  (un* ,  vn'  ,6n2 )  — •  ( u,v,d )  as  n  —  oo.  Recall  the  proof  of  Theorem  2.1  and  let 
Hn  denote  the  finite  dimensional  subspace  of  H  given  by  IIn  =  7/"1  X  H^2 ■  Note  that  II n  C  V  for 
all  n  =  (n1,n2)  6  Z+  x  Z+.  For  A  >  0,n  =  (ni,ri2)  G  Z+  X  Z+,  and  q  £  Q,  define  A"  (q)  £  C(IIn) 

by 

Anx(q)  = 


A2  +  AC711  (q)  +  A"'  (q)  Ln  (<?)* 

-Ln(q)  (l+A-M,Ba(?))J 


It  follows  that 


(A£  (?)  *",¥")  =  (AA(g)*B,tf") 


(3.12) 
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for  $n  =  ($” 1 ,  $£2)  i  ®  =  ($"* ,  $£2)  €  where  the  inner  product  (•,  •)  on  H  is  given  by  (2.13) 
and  the  operator  A\(q)  6  C(V,V*)  is  given  by  (2.15).  From  (2.16)  and  (3.12)  we  obtain  that 

</iX(9)*n,*n)  > /3a ||$12 ,  $ner  (3.13) 


*1“  =  P?'0  +  (Cni  (qn)  +  A)/V*V, 

*?2  =  A-1(P2n2T/+In(?n)P1nV), 

and  $  =  ($i,$2)  €  V*  is  given  by 

$1  =  V’  +  (C’(<?)  +  A)v>, 

$2  =  A-1  (??  +  L  (q)  ip) . 

Using  assumptions  (C),  (E),  and  (G),  it  is  not  difficult  to  argue  that  the  following  assertions 
hold. 

i) 

lim  ||$n  -  $||.  =  0, 

n  — *oo 


where  ]|-||.  denotes  the  usual  operator  norm  on  V*. 

(ii)  The  map  q  — ►  A\ (q)  $  is  continuous  from  Q  C  Q  into  V *  for  each  $  £  V,  and  in  particular 


lim  HAA(qn)V  -  A*(q)V ||.  =  0, 

n — *oo 


$  e  v 


whenever  qn  — ►  q. 


(iii)  There  exists  a  positive  constant  px,  independent  of  q  €  Q,  for  which 

||4a(?)*||.<pa||*||,  *  e  F. 

Furthermore,  if  we  let  Pn  :  H  — ►  Hn  denote  the  orthogonal  projection  of  H  onto  Hn  with  respect 
to  the  inner  product  (•,•)  given  by  (2.13),  then  =  (P”1  $1, P2n2*2) ,  *  =  (*i,*2)  €  P.  It 
follows  from  assumption  (G)  that 

lim  HP71#  -  ¥||  =  0,  *  e  V,  (3.16) 

n — ►  00 

and 

lim  IP71*  -  *1  =  0,  fell. 

n— ►  00 

Setting  =  ( u,6 )  and  =  (u"1,#"2),  from  (3.13)  and  (iii)  above,  we  obtain 
/?A||*n  -  Pn$||2  <  (4)((9n)(*n  -  P71*),*"  -  Pn'V) 

=  (Anx  ( qn )  Vn  -  Ax  (q)  *n  -  Pn*) 

+  (Ax(q)V  -  Ax(qn)*,*n  -  Pn*} 

+  (Ax  ( qn ) #  -  4.x  ( qn )  Pn«,  «n  -  P71* ) 

<  ll*n  -  *||.  II*"  -  P"*||  +  ||4a  (?)  *  -  ^  (?n)  *||. 
x  ||*"  _  p«*||  +  ||4a  (?n) (*  ~  PnV)H.  ||*n  -  Pn*|| 

<  ^  { ll*n  ~  *H2.  +  II^a  (?)  9 -A,  (?**)  *||a.  +  p\  ||*  -  Pn*||2 } 

+  yll*n  -  Pn*||2. 

It  follows  therefore,  from  (i),  (iii)  and  (3.16)  above,  that 

Urn  H*7*  -  P"*||2  <  i  lim  { ||$n  —  $||2  +  || 4 a  (9)  'P  —  4x  (?n)  *||2 

71-700  71-00  (3.17) 

4  p\  II*7*  —  Pn*||2  }  =  0. 
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The  estimate  (3.17)  together  with  (3.14),  (3.15),  (3.16)  and  assumption  (G)  yield  that 

lim  (uBl,i;n*,0Ba)  =  (u,w,0) 

n— »oo 

in  Ar  and  the  convergence  asserted  in  (3.9)  has  been  established. 

Remark  We  note  that  if  the  admissible  parameter  set  Q  is  infinite  dimensional  (i.e.  functional), 
its  discretization  could  also  be  included  within  our  general  approximation  framework  and  conver¬ 
gence  theory,  as  well.  For  a  detailed  description  of  how  this  can  be  accomplished,  see  for  example, 
[3], 


4.  Examples  and  Numerical  Results 

We  consider  the  basic  equations  of  one  dimensional  linear  thermoelasticity  (see  [5]  [16]) 

d2u  d3u  d2u 

PW  ~  Qv  (A  +  2/i)  dtfdt2  ~(X  +  2p}dg2 
06 

+  aT(3A  +  2p)—  = /0,  0  <  77  <  C,  t>0 

dO  d29  -  d3u 

pcm~KW+  Ti  +  2p)Wdi=90' 


(4.1) 


(4.2) 


0  <  T)  <  t,  t  >  0. 

With  the  introduction  into  the  first  equation  above  of  the  Voigt-Kelvin  viscoelastic  damping  term, 
equations  (4.1)  and  (4.2)  describe  the  longitudinal,  or  axial,  vibrations  of  a  thin  visco-thermoelastic 


rod  of  length  t.  In  the  above  equations,  u  denotes  the  axial  displacement,  6  the  absolute  tempera¬ 
ture,  p  is  the  mass  density  of  the  rod,  A  and  p  are  the  Lame  parameters,  c  is  the  rod’s  specific  heat 
and  k  is  its  thermal  conductivity.  The  positive  constant  9  is  known  as  the  reference  temperature  - 
the  absolute  temperature  of  a  stress  free  reference  state  for  the  rod.  The  functions  /o  and  go  repre¬ 
sent,  respectively,  an  externally  applied  axial  force  and  thermal  input.  The  nonnegative  constants 
ay  and  oj  denote  the  viscosity  coefficient  and  the  coefficient  of  thermal  expansion,  respectively. 
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We  are  interested  in  studying  the  partial  differential  equations  (4.1),  (4.2)  together  with  initial 
conditions  of  the  form 


U  (0,  Tj)  =  U0  (7?)  ,  —  (0,  7?)  =  U0  (7/)  ,  0  (0, 7?)  =  0O  (t?)  ,  (4.3) 

for  0  <  T]  <  £,  and  one  or  the  other  of  the  following  two  sets  of  boundary  conditions. 

(I)  Clamped  ends  with  temperature  fixed  - 

u  (f,0)  =  ii  (t,i)  =  0,  <>0 

0  (1,0)  =  0(t,C)  =  0. 

(II)  Clamped  and  insulated  ends  - 


u  (f,0)  =  u(t,£)  =  0,  t  >  b 

89  69 

^«,0)  =  ^(M)  =  0,  <>0. 

We  note  that  our  theory  can  of  course  handle  a  variety  of  other  boundary  conditions.  For  our 
purposes  here,  however,  to  simply  illustrate  the  application  of  our  general  appraoch,  the  two  given 
in  (I)  and  (II)  above  will  suffice. 

We  consider  the  problem  of  identifying  some  subset  of  the  parameters 

q  =  (av,Qr,A,/i,K)  (4.4) 

over  some  closed  and  bounded  subset  Q  contained  in  the  positive  orthant  of  the  Euclidean  space 
Q  =  R5.  Once  again  we  note  that  our  theory,  or  some  relatively  minor  modification  of  it,  can  handle 
inverse  problems  involving  the  estimation  of  any  of  the  parameters,  input  functions,  or  initial  data 
appearing  in  either  the  differential  equations  (4.1),  (4.2)  or  initial  conditions  (4.3).  We  consider 
only  the  subset  given  in  (4.4)  simply  for  the  purpose  of  illustration. 
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In  order  to  apply  the  abstract  theory  developed  in  the  second  and  third  sections  above,  we  make 
the  following  definitions,  identifications  and  assumptions.  We  set  H\  =  ^(O, t)  endowed  with  the 
standard  £2  inner  product  (to  be  denoted  by-{-,-)i)  and  set  H\  =  £2(0, f)  with  the  weighted  inner 
product 

e 

{'Pi'I’h  =  j 
0 

In  the  case  of  the  boundary  conditions  (I)  we  define  both  V\  and  Vo  to  be  the  Sobolev  space 
Hq(Q, l)  with  the  boundary  conditions  (II)  we  take  V\  =  Hq( 0,1)  and  V2  =  Hx(0,l).  The  Sobolev 
spaces  f/o(0,£)  and  are  assumed  to  be  endowed  with  the  usual  Sobolev  inner  products  and 

induced  Sobolev  norms. 

For  each  q  £  Q  we  define  the  operators  Aj(q)  £  C  ( Vj ,  V’ )  ,j  =  1,2,  C  (q)  £  £(Fi,Vj*),  and 
L(q)e£(  V2,V2')by 


1 

(A\  (q)<fi)(i>)  =  /  —  2fiDipDrp,  q>,il>eV 1, 

J  P 


(M  (q)v)W 


f. 

(C(q)<p)W  =  J 


t 

=  - 

J  pc 


DipDip,  <p,tl>€V2, 


av  iA  +  £)(pDifr,  ip,xp  £  Vi, 
P 


and 


(^{q)V>)W  =  ~  [  <*T - £VV>,  £  Vi,0  €  V2- 

J  pc 

0 

Definition  (2.4)  then  yields  that  the  operator  L(q)m  £  £(V 2,  Vj*)  is  given  by 

e 

{L(q)M  ip)  (<fi)  =  J  Dip<p,  <peVi,ip€V2. 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


We  assume  that  u0  £  Hq(0,£),  vo,0q  £  L2(0,£),  fo,go  €  £1  (0, 7,  L2(Q,£))  and  set  /  (t,rj)  = 
-pfa  (f,  77),  g  (t,T])  =  -Lg0  (<,7?)  for  almost  every  (t,  q)  £  [0,T]  X  [0,£].  It  is  then  a  simple  and  straight 
forward  matter  to  show  that  the  assumptions  (A)  -  (F)  stated  in  Section  2  are  satisfied. 
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For  each  N  =  1,2,...  let  }^_0  denote  the  standard  linear  5-spline  (i.e.  “pup- tent”  or  “hat”) 
functions  on  the  interval  [0, £]  defined  with  respect  to  the  uniform  mesh,  {Q,l/N,2£/N, . . .  ,^}. 
If  for  Tii  =2,3,...  we  set  5"1  =  span  {yjJ1  1  and  for  each  ri2  =  2,3,...  we  set  H £2  - 
span  1  in  the  case  of  boundary  conditions  (I)  or  H^2  =  span  {pj2  }n^.0  in  the  case  of 

boundary  conditions  (II),  then  H J’  C  V},  j  =  1,2,  and  ,  using  standard  approximation  estimates 
for  interpolatory  splines  (see  [14])  together  with  the  Schmidt  inequality  (see  [11])  it  is  not  difficult 
to  argue  that  assumption  (G)  is  satisfied  (see,  for  example,  [10],  [11]). 

With  the  spline  bases  for  H "*  and  H^2  defined  in  the  previous  paragraph  and  using  the  definitions 
(4.5)  -  (4.9),  it  is  straight  forward  to  write  the  finite  dimensional  initial  value  problems  (3.1)  -  (3.3) 
(or  equivalently,  (3.4),  (3.5))  in  matrix  form.  In  the  examples  to  follow  we  took  the  observation 
space  Z  to  be  Z,2(0,^)  and  chose  the  functionals  as 

e 

$i(x,q;z)  =  ((«,  v,  6) ,  q\  z)  = 

o 

i  =  1,2,...  ,v. 

We  used  the  IMSL  routine  ZXSSQ ,  an  implementation  of  the  iterative  Newton’s  method /Steepest 
descent  hybrid  Levenberg-Marquardt  algorithm  to  solve  the  resulting  finite  dimensional  nonlinear 
least  squares  minimization  problems  (IDn).  The  IMSL  routine  IVPAG,  an  implementation  of  Gear’s 
method  for  stiff  systems,  was  used  in  each  iteration  to  solve  the  initial  value  problems  (3.4),  (3.5) 
for  a  given  choice  of  q  6  Q,  whenever  required.  A  composite  two-point  Gauss  Legendre  quadrature 
formula  was  used  to  numerically  compute  integrals  whenever  necessary  (in  inner  products  when 
computing  orthogonal  projections,  for  example).  All  codes  were  written  in  Fortran  and  all  com¬ 
putations  were  carried  out  on  the  Cray  XMP/48  at  the  San  Diego  Supercomputer  Center.  We 
note  however,  that  since  we  were  only  interested  in  estimating  constant  (as  opposed  to  functional) 
parameters,  the  use  of  the  supercomputer  was  probably  superfluous. 


19 


Example  4.1  In  this  example  we  seek  to  1)  simply  illustrate  that  our  general  approach  is  valid 
in  practice  as  well  as  in  theory  and  that  the  spline  based  scheme  that  has  been  proposed  above 
performs  satisfactorily,  and  2)  to  obtain  some  feel  for  how  good  of  a  fit  can  actually  be  obtained. 
This  second  goal  will  aid  us  to  a  certain  extent  in  interpreting  our  findings  in  Example  4.2  to  follow 
where  we  consider  a  somewhat  more  realistic  inverse  problem. 

In  order  to  achieve  the  two  above  stated  goals  the  following  procedure  was  used.  So  called  “true” 
values  for  the  unknown  parameters  and  functions  u  and  6  were  chosen.  Then  /0  and  go  were 
determined  so  that  the  partial  differential  equations  (4.1)  and  (4.2)  were  satisfied  and  Uo,uo,  and 
60  were  chosen  in  accordance  with  (4.3).  In  this  example  we  considered  the  boundary  conditions 
specified  in  (I). 

We  chose 

p  =  ay  =  X  =  p  =  aj  =  c  =  n  =  6  =  i  =  l, 
u(t,T])  =  (sin  t)  (sin  7r7?) ,  0  <  t]  <  1,  t  >  0, 

and 

6(t,ri)  =  el  (ri2  -  T])  ,  0  <  rj  <  1,  t  >  0. 

Observations  were  generated  according  to 

z(ti,i))  =  (sinf,-)(sin7T7?),  0  <  rj  <  1, 

for  ti  =  .5 i,i  =  1,2,. . .  ,8  =  v.  The  discretization  levels  n\  and  in  our  spline  scheme  were  always 
chosen  the  same.  Thus,  without  fear  of  confusion,  we  shall  refer  to  their  value  as  n  =  n\  =  n?. 
Note  that  in  this  case,  the  dimension  of  the  approximating  system  of  ordinary  differential  equations 
(3.4),  (3.5)  is  3n  -  3 

(a)  In  this  example  we  seek  to  identify  q  =  (qv,pt)-  As  an  initial  guess  with  which  to  start  the 
iterative  optimizing  search  we  used  q°  -  (.25,. 25).  Our  results  are  summarized  in  Table  4.1a 
below. 
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n 

r*n 

Oty 

Qj 

Jn{q°) 

Jn(qn) 

CPU(m:s) 

4 

.99045 

.62424 

3.41 

5. 13x10-' 

0:6.34 

8 

1.0159 

.90622 

3.73 

mmsum. 

1:46.17 

12 

1.0025 

.98411 

3.88 

KBfiSjiiil 

5:46.50 

16 

1.0020 

1.0037 

3.94 

11:45.45 

I 


Table  4.1a 


(b)  In  this  example  we  consider  a  somewhat  more  challenging  problem;  we  seek  to  identify  the  five 
parameters  q  —  (qv/,  q^,  A,/x,  k).  As  an  initial  guess  we  used  q°  =  (.25,  .25,  .25,  .25,  .25).  Our 
results  are  displayed  in  Table  4.1b  below 


n 

Qj 

An 

A" 

kn 

Jn(q°) 

Jn(qn) 

CPU(m:s) 

m 

.89827 

1.0720 

.99102 

1.0393 

.93712 

8.25xl02 

2.24x10-' 

1:18.14 

8 

1.0090 

1.0102 

.94168 

.98898 

.98229 

4:56.44 

12 

.99566 

.99349 

.98076 

1.0016 

.99094 

■iltEEW 

27:42.41 

16 

1.0028 

1.0072 

.97233 

1.0058 

.99443 

40:19.79 

Table  4.1b 

Example  4.2  In  this  series  of  tests  we  indicate  how  our  approach  might  be  applied  in  practice, 
and  how  it  might  perform.  In  particular  we  illustrate  how  it  may  serve  as  an  aid  in  the  identification 
of  internal  dissipation  (i.e.  damping)  mechanisms  in  flexible  structures. 

We  consider  a  long,  slender  aluminum  (AL  T2024  -  T4)  rod  of  length  100  in.  which  is  clamped 
and  insulated  at  both  ends  (i.e.  boundary  conditions  (II)).  We  assume  that  it  is  set  to  vibrating 
via  a  mechanical  input  of  the  form 


fo{t,V) 


f  aoX[eue2]  (i?)sin37rt,  0  <  77  <  £,  0  <t<T0 

1  0,  To  <  t  <  T, 


where  \i  denotes  the  characteristic  function  for  the  interval  /  C  R.  We  assume  that  there  is  no 
externally  applied  thermal  input  -  i.e.  we  take 


9o{t,r])  =  0,  0  <  Tj  <  l,t  >  0, 
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and  we  assume  that  the  system  is  initially  at  rest.  That  is,  we  take  uo  =  t>o  =  0,  and  <?o  =  0,0  < 
7/  <  l. 

After  appropriate  normalization,  the  laboratory  or  experimentaJly  measured  values  for  the  ma¬ 
terial  parameters  appearing  in  the  system  (4.1),  (4.2)  are  given  by 

e  =  i 

p  =  9.82  x  10-2  aj  =  1 .29  x  10~3 

A  =  2.0635  x  10"1  C  =  5.39611  x  10"1  (  1U) 

p  =  1.1111  x  10"1  k  =  7.02176  x  10“7. 

There  is  no  generally  accepted  means  to  explicity  measure  the  viscosity  coefficient  o\-  in  the  labo¬ 
ratory.  Typically  its  value  is  determined  via  an  identification  procedure  such  as  the  one  which  we 
have  developed  here.  We  took  0  =  6.8  x  10*. 

With  the  boundary  conditions  given  in  (II),  the  natural  mode  shapes  for  the  homogeneous  sys¬ 
tem  corresponding  to  (4.1),  (4.2)  are  of  the  form  $o  =  (uo,0o)  =  (0,1),  =  (um,flra,j)  - 

(sin  cos  ^j^),  m  =  1,2,...  J  =  1,2,3,  where  for  each  m  -  1,2,...,  and  j  =  1.2,3, 

Cm j  is  a  complex  constant.  Thus  when  a  Galerkin  scheme  using  A/  +  1  elements  of  the  form 

=  (sin  cos  ,m  =  0,1,2 _ ,M,  is  employed  to  discretize  (4.1),  (4.2),  the  resulting 

finite  dimensional  equivalent  first  order  system  consists  of  M  -f  1  decoupled  first  order  systems,  one 
of  dimension  one  and  the  remaining  M  each  of  dimension  three.  We  used  such  a  scheme  with  the 
true  values  of  the  parameters  and  M  =  24  to  generate  simulated  observational  data  upon  which 
to  base  our  fits.  Once  again  in  our  spline  scheme  we  use  n  =  n\  -  n-i  to  denote  the  level  of 
discretization. 

The  presence  of  either,  and  certainly  both,  the  Voigt-Kelvin  viscoelastic  term  (i.e.  ov  >  0) 
and  the  thermoelastic  term  (i.e.  aj  >  0)  yield  some  form  of  mechanical  energy  dissipation.  In 
fact,  when  the  boundary  conditions  are  as  given  in  (II),  with  a  simple  change  of  variable  and  the 
construction  of  an  appropriate  Lyapunov  functional  (or,  equivalent  norm),  it  can  be  shown  that 
either  the  Voigt-Kelvin  damping  or  the  thermoelastic  coupling  leads  to  a  uniformly  exponentially 
stable  open-loop  or  unforced  system  (i.e.  |5(*;?)|  <  M  exp(u>t)  ,t  >  0,  with  u>  <  0,  where 
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:  t  >  o|  is  the  semigroup  corresponding  to  the  transformed  system  (see  [6]  and  [7]).  The 
simple  change  of  variable  consists  of  transforming  the  absolute  temperature  6  to  the  temperature 
variation  from  the  time  invariant  average  temperature  over  the  length  of  the  rod,  dAv  e-  That  such 
a  constant  average  temperature  exists,  independent  of  the  input  /o  and  initial  conditions  (4.3),  is  a 
direct  consequence  of  the  insulated  ends.  Indeed,  integrating  equation  (4.2)  from  0  to  f,  integrating 
by  parts  and  then  imposing  the  boundary  conditions  (II),  yields  the  desired  conclusion.  We  note 
that  even  without  the  change  of  variable,  the  mechanical  vibrations  of  the  unforced  system  still 
tend  to  zero  exponentialy  fast.  (In  this  case  the  absolute  temperature  tends  to  9ave  exponential)’ 
fast.)  Thus  the  presence  of  either  the  viscoelastic  (Voigt- Kelvin)  or  the  thermoelastic  effects  will 
exponentially  damp  unforced  mechanical  vibrations.  In  the  series  of  examples  to  follow  we  illustrate 
how  our  scheme  and  general  approach  might  be  used  to  aid  in  the  identification  of  the  mechanism 
for,  and/or  an  appropriate  model  for  an  observed  dissipation  of  mechanical  energy  in  a  flexible 
structure. 


(a)  In  this  test  we  took  a0  =  l,fi  =  .4,  s2  =  -435,  To  =  .5  and  T  =  5.  We  set  or  =  0  (i.e. 

no  Voigt-Kelvin  damping),  generated  observations  at  times  t,  =  .25 1,  i  =  0,1,2 . 20  via  the 

modal  scheme  described  above,  and  used  our  scheme  to  see  if  we  could  estimate  or-  All  other 
parameters  were  taken  as  given  in  (4.10)  and  c*v  was  held  fixed  at  zero.  For  an  initial  guess 
we  set  a®  =  10-2.  Our  results  are  summarized  in  Table  4.2a  below.  It  is  clear  from  the  table 
that  n  had  to  be  taken  quite  large  in  order  to  obtain  an  accurate  estimate  for  07-.  It  is  likely 
that  the  reason  for  this  is  that  thermoelastic  effects  in  metals  are  relatively  slight.  Consequently 
we  would  expect  that  the  thermo-mechanical  coupling  as  described  by  the  parameter  a 7-  is  not 
readily  discernable  by  our  approximation  and  identification  procedures. 
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n 

Gj- 

Jn(q°) 

Jn(qn) 

CPU(m:s) 

4 

3.0409xl0-3 

1.189x10-' 

1.187x10-' 

0:4.94 

8 

BSEliSS&BH 

1:11.12 

12 

esmbbh 

1  SEBI 

3:15.86 

16 

ffiillllSH 

HKgjSBIl 

21:32.26 

20 

4.471xl0-2 

45:36.64 

Table  4.2a 


(b)  This  test  is  similar  to  the  previous  one  described  in  (a)  above,  however,  this  time  we  introduced 
some  Voigt-Kelvin  damping  into  the  system.  We  set  ay  —  10-2  when  the  observations  were 
generated.  This  time  we  also  took  «o  =  5.0, T  -  3.0,  and  once  again  set  q“  =  10-2.  Our 
findings  are  reported  in  Table  4.2b.  In  general,  the  presence  of  the  Voigt-Kelvin  damping  tended 
to  make  our  task  somewhat  simpler  and  to  improve  the  overall  performance  of  our  scheme.  This 
is  not  surprising  since  Voigt-Kelvin  is  a  rather  strong  form  of  dissipation  and  its  presence  would 
be  likely  to  significantly  increase  the  stability  and  accuracy  of  our  numerical  integrator. 


n 

a'} 

Jn(qn ) 

CPU(m:s) 

B 

2.1998xl0-4 

1.607x10-' 

1.570x10-' 

0:54.21 

8 

BlfTfini'MI 

1:3.82 

16 

bessebh 

■OTOlRIiBM 

4:35.93 

20 

wimnrBB 

14:17.86 

Table  4.2b 


(c)  In  this  example  we  first  attempt  to  fit  a  model  with  only  thermoelastic  damping  to  observations 
of  a  system  which  has  both  visco-  and  theoremoelastic  dissipation.  We  then  try  to  fit  a  model 
with  viscoelastic  damping  only  to  a  system  in  which  the  dissipation  is  due  only  to  thermal  effects. 
In  both  tests  the  input  excitation  and  observations  are  as  they  were  in  test  (b)  above. 

For  the  first  test  we  generated  observations  with  av  =  10-2  and  aj  as  given  in  (4.10).  We 
then  attempted  to  estimate  the  parameter  oV'  with  aj  held  fixed  at  zero.  The  initial  guess  for 
was  taken  to  be  =  2.5  x  10_‘.  Our  results  are  given  in  Table  4.2c. 1. 

In  the  second  test,  observations  were  generated  with  -  0  and  a /  as  given  in  (4.10).  We 
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then  again  attempted  to  estimate  ay  with  ar  fixed  at  zero.  The  initial  guess  was  a°v  =  10  2 . 
The  results  for  this  test  are  given  in  Table  4/2c.2 


n 

Oty 

Jn(Q°) 

Jn(qn) 

CPU(m:s) 

n 

7.1179xl0~3 

4.234xl0-1 

1. 466X10"1 

0:10.23 

8 

1.2178x10-" 

4.122x10-* 

4.863x10“" 

0:44.62 

16 

9.7023xl0"3 

4.157x10-* 

1.064x10"" 

9:12.97 

20 

9.7812X10"3 

4.165x10-* 

1.135x10-" 

16:17.10 

Table  4.2c.l 


n 

*nv 

Jn(Q°) 

Jn(qn) 

CPU(m:s) 

4 

1.4501xl0"3 

3.329xl0-1 

3.019X10"1 

0:8.68 

8 

IfMWiffliEl  liaM 

‘Btffi'irn'Mi 

0:24.14 

16 

IKSII 

IFTT5TgHB 

3:43.28 

20 

Tyff'TffnvfiM 

12:53.48 

Table  4.2c.2 

,  For  the  system  that  we  have  considered  here,  the  results  given  in  Tables  4.2c. 1  and  4.2c.2  seem 

|  to  contribute  to  the  conclusion  that  the  effects  of  thermoelastic  damping  are  insignificant.  Indeed, 

in  the  first  test,  a  relatively  accurate  estimate  for  the  viscosity  coefficient  ay  can  be  obtained  with 

I  a  model  that  ignores  all  thermoelastic  effects.  The  second  test,  although  not  conclusively,  appeared 

■ 

1  to  indicate  that  the  presence  of  thermoelastic  effects  in  the  data  will  have  only  a  minimal  effect  on 

j  the  estimates  obtained  for  the  viscosity  term  ay 

!  We  note  at  this  point  however,  with  a  word  of  caution,  that  the  interpretation  given  in  the  pre¬ 

vious  paragraph  is  not  intended  as  a  sweeping  generalization  and  conclusion  about  thermoelastic 
effects  in  the  vibration  of  flexible  structures.  Indeed,  we  present  these  findings  simply  as  an  illus¬ 
tration  of  how  our  general  approach  and  schemes  could  be  used  in  such  a  damping  study.  To  draw 
such  conclusions  would  necessarily  require  an  exhaustive  investigation  involving  both  simulation 
|  and  computational  studies  as  well  as  laboratory  experimentation. 
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5.  Summary  and  Concluding  Remarks 

In  this  paper  we  have  considered  inverse  or  parameter  identification  problems  for  thermoelastic 
systems.  We  defined  a  class  of  abstract  infinite  dimensional  systems  consisting  of  a  second  order 
elasticity  equation  coupled  with  a  first  order  abstract  parabolic  equation.  We  established  the 
well-posedness  of  this  system  by  rewriting  it  as  an  equivalent  abstract  first  order  system  set  in  an 
appropriate  product  Hilbert  space  and  then  applying  results  from  linear  semigroup  theory.  A  class  of 
identification  problems  in  which  the  state  constraints  were  given  by  the  mild  solution  to  the  abstract 
thermoelastic  system  were  defined.  We  then  proceeded  to  develop  an  abstract  approximation 
framework  based  upon  generic  Galerkin  approximation  of  the  mild  solutions.  Using  approximation 
results  from  linear  semigroup  theory,  we  were,  under  appropriate  assumptions,  able  to  establish  that 
the  solutions  to  the  resulting  finite  dimensioned  approximating  parameter  estimation  problems,  in 
some  sense,  approximate  a  solution  to  the  original  infinte  dimensional  identification  problem.  We 
demonstrated  that  our  approach  could  be  applied  to  the  basic  equations  of  linear  thermoelasticity 
and  provided  computational  results  as  evidence  that  our  schemes  are  not  only  theoretically  well 
founded,  but  numerically  sound  and  pactical  as  well. 

Many  interesting  and  important  open  questions  related  to  the  identification  and  control  of  ab¬ 
stract  thermoelastic  systems  remain.  In  particular,  nonlinear  systems  and  especially  ones  in  which 
the  material  parameters  are  permitted  to  be  temperature  dependent  need  to  be  considered.  Also, 
numerical  studies  using  actual  experimental  data  rather  than  simulation  data  should  be  carried 
out  and  reported  on.  With  regard  to  the  optimal  control  problem,  a  complete  understanding  of 
the  asymptotic  behavior  of  the  spectrum  of  the  operator  A{q)  given  by  (2.6)  would  be  extremely 
useful  as  would  it  be  for  just  the  basic  equations  of  linear  themoelasticity  (4.1),  (4.2).  Indeed, 
for  the  boundary  equations  (II)  it  can  be  shown  that  the  eigenvalues  asymptotically  approach  a 
straight  line  in  the  left  half  plane.  For  other  standard  sets  of  boundary  conditions,  in  particular 
the  Dirichlet/Dirichlet  case  given  in  (I),  this  author  knows  of  no  similar  types  of  theoretical  results. 
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Related  to  this  is  the  question  of  uniform  exponential  stability  of  the  open-loop  semigroups.  For 
some  sets  of  boundary  conditions,  this  can  be  established  (see  [7]).  However,  the  question  in  general 
is  still  open.  Finally,  when  doing  approximation  for  LQG  control,  among  the  conditions  required 
to  establish  convergence  of  the  optimal  closed-loop  feedback  gains,  is  the  preservation  of  uniform 
exponential  stability  under  approximation  (see,  for  example,  [6]).  Even  for  a  linear  spline  based 
Galerkin  scheme  for  the  basic  equations  of  linear  thermoelasticity  with  the  boundary  conditions 
given  in  (II),  this  condition  does  not  appear  to  hold.  Schemes  for  which  this  condition  can  be 
verified  (other  than  modal  schemes  of  course,  for  which  it  is  trivially  true)  need  to  be  developed 
and  studied. 
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